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I. THE GEIGER-MARSDEN EXPERIMENTS 


In 1911 Ernest Rutherford proposed a nuclear theory of the atom as 
a result of the alpha particle scattering experiments performed by Geiger 
and Marsden. In Rutherford's own words: , 


"In the early days I had observed the scattering of alpha 
particles, and Dr. Geiger in my laboratory had examined it 
in detail. He found, in thin pieces of heavy metal that 
the scattering was usually small, of the order of one de- 
gree. One day Geiger came to me and said, 'Don't you think 
that young Marsden, whom I am training in radioactive 
methods, ought to begin a small research?' Now I had 
thought that too, so I said, “Why not let him see if any 
alpha particles can be scattered through a large angle?' 

I may tell you in confidence that I did not believe that 
they would be, since we knew that the alpha particle was 

a very fast, massive particle, with a great deal of energy, 
and you could show that if the scattering was due to the 
accumulated effect of a number of small scatterings the 
chance of an alpha particle's being scattered backwards 
was very small. Then I remember two or three days later 
Geiger coming to me in great excitement and saying, ‘We 
have been able to get some of the alpha particles coming 
backwards...’ It was quite the most incredible event 

that has ever happened to me in my life. It was almost 

as incredible as if you fired a 15-inch shell at a piece 
of tissue paper and it came back and hit you. On consid- 
eration, I realized that this scattering backwards must be 
the result of a single collision, and when I made calcula- 
tions I saw that it was impossible to get anything of 

that order of magnitude unless you took a system in which 
the greater part of the mass of the atom was concentrated 
in a minute nucleus. It was then that I had the idea of 
an atom with a minute massive center carrying a charge. 

I worked out mathematically what laws the scattering 
should obey, and I found that the number of particles 
scattered through a given angle should be proportional to 
the thickness of the scattering foil, the square of the 
nuclear charge, and inversely proportional to the fourth 
power of the velocity. These deductions were later ver- 
ified by Geiger and Marsden in a series of beautiful 
experiments.''* | 


The experiments performed by Geiger and Marsden made use of properties 
of alpha particles previously discovered by Rutherford. Figure 1 is a 
sketch of Geiger's scattering apparatus. The alpha particles emitted by 
the radium in the cone-shaped glass tube at the left strike the metal foil 
placed over the slit. When the deflected alpha particles strike the zinc 


*E. Rutherford, "The Development of the Theory of Atomic Structure", from 
Background to Modern Science, New York: MacMillan & Co., 1940. 
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FIGURE | 


sulphide and glass scintillation screen at the right end of the tube, the 
flashes of light can be viewed through the telescope and the number of 
flashes recorded. 





_ Figure 2 is a sketch of the Geiger-Marsden experiment to investigate 
"diffuse reflection", or large-angle deflections of alpha particles. 


Alpha particles from a variety of sources were aimed at foils of lead, 
gold, platinum, tin, Silver, copper, iron, and aluminum. The scintillation 
Screen was placed at a 90° angle to the foil, and the lead block shielded 
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the screen from alpha particles coming directly from the cone. 
angle deflections were observed as before by means of a microsco 
behind the screen. 
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In these series of experiments Geiger and Marsden found that less 


than one alpha particle in 8000 was deflected back toward the sc 
the number of alpha particles scattered onto the screen increase 
increasing nuclear charge as well as increasing thickness of the 
that the more energetic the alpha particles, the higher the prop 
them scattered through a large angle. 
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Ii. THE PHYSICS BEHIND THE SCATER PROGRAMS 


A. SCATR1: THE EXPERIMENT | 


SCATRI calculates the number of 5.3 MeV alpha particles scattered at 
angles between 10° and 180° by using the Rutherford formula for the frac- 
tion of beam particles scattered: 


ee N D7 A 
D 16 L2 Sin*(@/2) 


where N is the number of target atoms in a square piece of foil lem on a 
Side; A is the area of one of the counters that detects scattered alpha 
particles (lem? here); L is the distance from the target to the counters 
(10cm here); 9 is the scattering angle; and D is the "distance of closest 
approach," a convenient expression that is given by: 


Z e2 
, Z, e 


D = E 


Here Z, is the charge of the alpha particle (2 electron charges); Z, is the 
charge of the gold target nucleii (79 electron charges); e is the charge on 
an electron; and E is the kinetic energy of the incident alpha particle 
(5.3 MeV here). 


N, the number of target atoms per square cm is calculated using 


V-H fT 
W 


where V is Avogadros number (6.02 x 1023 atoms per mole) ; H is the density 


> 
of the target foil (here 19.3 grams/cm%); T is the thickness of the foil 
(here 107-*cm or 1 micrometer) ;and W is the atomic weight of the target 
material (here 197 grams/mole). 
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The number of scattered alpha particles going into each counter is 
calculated by multiplying the calculated value of Py at the appropriate 9 
by the number of incident alpha particles: 


where B is the total number of alpha particles hitting the target during 
the counting period. Since 10° particles hit the target each second and 
the counting period is 10° second, B is 10/9 particles. 


The randomness of any scattering experiment is simulated by adding a 
normally distributed random variable with mean of zero and standard devia- 
tion of (N_) 2 to the value of N.. This simulates the statistics of 


actual particle detection. The printed values of the number of particles 
detected at each angle will show very little variation when the number of 
detected particles is large, but considerable variation when the number of 


detected particles is small. For large scattering angles at which the 
number of detected particles is only about 100, one standard deviation is (100) 1/2 


or 10 and successive measurements taken at the same angle will be quite dif- 
ferent. (The variation actually follows a Poisson distribution, but for 
numbers greater than about 20, the normal or Gaussian distribution is an 
often-used and accurate approximation.) 


The fraction of scattered particles printed for each angle is calcu- 
lated by dividing the number of scattered particles by the number of 
incident particles (19!9), 





B. SCATR2: THEORETICAL PREDICTION FROM THREE MODELS 


The Hard-Sphere Model 


The fraction of the incident particles scattered into a counter of area A 


at a distance L from the target by a target made up of atoms which are hard 
spheres of radius R is: 


Py R? NA 
D Te 


where p, is the probability of interaction between the incident alpha particles 
and one of the target spheres, (here po = 107°), and N is the number of spheres 
in a piece of target foil lcm on a side. The fraction scattered into each 
counter is completely independent of the scattering angle. Po is fixed at 10° 
and R is fixed at 107° em, the actual size of a typical atom (gold atoms are 
slightly larger, hydrogen atoms are slightly smaller). Since this model is 
wrong, there is no way of "knowing" what the real value of p, should be. If the 
probability of interaction is taken as 1, all of the incident particles will 
scatter from a target foil only one layer of atoms (107%em) thick. Py was 
chosen to give convenient numbers of scattered particles for the foil thickness 
used in the experiment simulation in SCATR1. Other values of Po would merely 
change the total number of scattered particles but not the angular distribution. 


As before, N is calculated by 


where V is Avogadro's number (6.02 x 107°); H is the target foil density: T is 
the target foil thickness; and W is the target foil atomic weight. Values of 
H and W for copper, silver, and gold are listed in the nomenclature table. 


The Thomson Model 


The fraction of the incident particles scattered into a counter of area A 
at a distance L from the target by atoms that are diffuse spheres of positive 
charge with electrons suspended within (Thomson atoms) is zero for scattering 
angles greater than a maximum angle 9. The maximum angle at which a Thomson 
atom can scatter an incident alpha particle of kinetic energy E is 


: ZZ e@ 
9 = eS ee | -1(D 


The combination Of conditions for which 9 is a maximum within the constraints 
of the program is: Z,=79 (gold), E=2MeV. Even for these values, Qo. is only 
0.03°, and no particles are scattered into any of the counters. 
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The Rutherford Model 


The fraction of incident particles scattered into a counter of area A at 
a distance L and angle 9 from a target made up of atoms in which a central 
charged nucleus creates a Coulomb potential (the Rutherford Model) is 


ND* A 


rk. = ——_ 
D 16 = sin'(@/2) 


where N is the number of target atoms per square centimeter, calculated by 


e 
VHT 
, RE i Bea ea ° 
° W 
Again, V is Avogadro's number and H, W, and T are the target foil density, 
atomic number, and thickness. 
é D is given, as before, by 
Zz 2. e* 
os oe ee 
E 
Z, and E are the charge and kinetic energy of the alpha particle, Z, is the 
charge of each target atom nucleus, and e is the electron charge. 
The terrible problem of units is avoided here and in most of these 
calculations by replacing e* by the famous expression 
e* =r mc? - 
ee 
r, is the "classical electron radius": 
: ay Fa" 
‘ * 210-3230 meters = 2.8 femtometers 
and m,c*, the rest energy of an electron, is 
= - 


ac* = ~511 MeV 





The shape of the angular distribution is determined by sin'(9/2). The 
number of particles scattered is proportional to the thickness of the Foil, 
the area of each counter, the density of the foil, and the square of the atomic 
number of the foil atoms, The number of particles scattered is inversely 
proportional to the atomic weight of the target atoms, to the square of the 
alpha particle energy and to the square of the distance from the foil to the 
counters. 


When § approaches zero, the fraction of particles scattered at the angle 
9 gets very large. The fraction is prevented from going to infinity (it, in 
fact, stays less than 1) since for very small 9, the range of 8 measured by a 
counter must be very small and hence A, the area of the counter , must be very 
small. Although 1/sin*t(@/2) goes to infinity for 9=0, the area of the 
counter at 9=0 is also zero, and N_ remains finite. For small angles, where 
the fraction scattered is relatively high, each alpha particle has a fairly 
high probability of being scattered by a second atom on its way out of the 
foil, and the scattering data at angles less than a few degrees (not meas- 
ured in this simulation) will reflect this multiple scattering. 




















C. SCATR3S: TRAJECTORIES 


The trajectory of an alpha particle being scattered by a stationary 
charged nucleus is described (in polar coordinates: r, 6) by 


oe D (See, for instance, 
Xo +350r = Cos92* 2 2 sin = 0. E. J. Konipinski, 
CLASSICAL DESCRIPTIONS 
OF MOTION, p.98.) 


Here x is the perpendicular distance between an extension of the alpha 
particle's initial path and the atomic nucleus (the ‘impact parameter'), 

D, the "distance of closest approach,' is how close the alpha particle could 
come to the nucleus if its initial path were aiming straight for the nucleus. 
D is given by 


zZZ,e° 
}y ae eee. 
E 


where Z, and E are the charge and kinetic energy of the alpha particle, Zo 

is the charge of the nucleus, and e is the charge on an electron. 
Converted to rectangular coordinates, with the alpha particle arriving 

from the negative y direction, the trajectory equation is 

i 

Rie 


D ‘ 
x2 + —(x2 4. y*) 5 


3 34 


y~- Xx ’; 0 


where 


x = r sind : y = -r cos® 


When solved for x as a function of y, this equation gives 


—_— 2: ae ee z 
ge x - @)* tg VP t= Pay 





The behavior of x as a function of y depends on the value of x_, the perpen- 
dicular distance between the initial alpha particle path and the nucleus: 


a) When Xo is greater than D/2: | _ : 
Only the positive sign in the equation for x is used: 


x xD i 


oe eo gare 
"2. OF = + 5 w+ |y a. Dy]* ) 
Oo 2 





The scattering angle is less than 90° and each value of y on 
the trajectory corresponds to only one value of x. 





b) When x, is less than D/2: 


Use of both the positive and negative sign in the expression 

for x gives two different values of x for each value of y. 

This is necessary, since the scattering angle is now greater 
than 90° and the outgoing alpha particle trajectory again has 
negative values of y. In this case, the maximum (least negative) 
value that y attains is 


i 
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c) When x, equals D/2: 


The trajectory equation becomes 
1 


at (x? + y2)? + y- x = 0 
and x as a function of y is 
: D/ 2y ‘2 
2\ 2y +D a. 
Now the maximum value that y takes on is 
7 


In each case 9, the angle through which the alpha particle is 
scattered, is given by 


tan(@/2) = 
2 


xX st 
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III. PROGRAM LISTINGS (Seatrl, Seatr2, Scatr3) 


REM SCATR1 - ALPHA PARTICLE SCATTERING EXPERIMENT SIMULATION 

REM COPYRIGHT 1973 - STATE UNIVERSITY OF NEW YOR 

REM DEVELOPED BY Je BRRESMANs Ae CAGGIANOs 4ND De SCARL 

REM PROGRAMMED BY Ce LOSIK AND De SCARLSAREVISED RY Me WEISNER 
REM SCATER IS 4 SET OF ASSOCIATED PROGRAMS CONTAINING SCATRI; 
REM SCATR@s, AND SCATRS 

REM LATEST REVISION: 12-12-73 

RANDOMIZE 

LET NAA 

DIM GC5)sHC3)eWVOG)92463) 

REM AVOGADRO'S NUMBER 

LET V4=6-92EF23 

REM ELECTRON REST ENERGY (MEV) 

LET B4=.511 

REM CLASSICAL ELECTRON RADIUS (CM) 

LET D4=2.e8k-13 

REM 4SLPHA PARTICLE CHARGE 

LET.4&4=2 

REM ATOMIC NUMBER QF GOLD 

LET 2¢€3)=79 

REM ATOMIC WEIGHT OF GOLD 

LET WC¢C3)=197 

REM DENSITY OF GOLD 

LET HC3)=19.-3 

REM LOG OF NO. OF REAM PARTICLES 

BES ees oe 

LET R=1A4 

Lt fF AS Se 3 

LET T=1 

LET Tas lE<5*T 

LET AA4=1 

LET L4=14 

PRINT “INSTRUCTIONS (93=N0s1=YES)"$ 

INPUT If | 

IF [=0 THEN 580 

IF I<>1 THEN 420 | | 

PRINT'A VERY RADIOACTIVE POLONIUM SOURCE PRODUCES 4 REAM OF" 
PRINT'19%,009% ALPHA PARTICLES PER SECONDs BACH HAVING AN ENERGY 
PRINT™'OF 5-3 MEVe THE REAM STRIKES A 1 MICROMETER THICK GOLD FOIL." 
PRINT"FIVE ALPHA PARTICLE DETECTORS CAN RF SET UP TO COUNT THE" 
PRINT'NUMBER OF SCATTERED PARTICLES DURING 4 COUNTING PERTOD OF 
PRINT™“199,009 SECONDS CAROUT 28 ROURS)Ie EACH DETECTOR HAS 4N ARFA™ 
PRINT™OF 1 SQ CM AND IS 19 CM SWAY FROM THE GOLD FOIL" 
PRINT'THE ANGULAR LOCATION OF FACH OF THE COUNTERS CAN RE SET" 
PRINT'RETYEEN THE LIMITS OF 14 DEGREES AND 189% DEGRESS RY TYPING" 
PRINT™A NUMRER BETWEEN 19 SND 18% WHEN THE COMPUTER 4S5”S: CC)=." 
PRINT’ZERO DEGRERS IS FORWARD SCATTERING AND 1890 DEGREES 185 " 
PRINT'RACKWOARD SCATTERINGe TYPE CO)=4 TO FND THR PROGRAM." 


30 REM TARGET PARTICLES PER Cro 


LET N=T4*V4*HC3)7WC3) 
REM DISTANCE OF CLOSEST APPROACH 
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SCATR2 


19% REM SCATR® - ALPHA PARTICLE SCATTERING - THFORFTICAL PREDICTIONS 
114 REM COPYRIGHT 1973 - STATE UNIVFRSITY OF NEW YOoRX ; 
{29 REM MODELS DEVELOPED BY Je BRESMANs Ae CAGGIANOs 4ND De SCARL 
1.2% REM PROGRAMMED RY Ce LOSIK AND De SCARL» REVISED RY Me WEISNER 
14% REM SCATER IS 4 SET OF ASSOCIATED PROGRAMS CONTAINING SCATRI> 
159 REM SCATR2s AND SCATR3 
e 16% REM LATEST REVISION: 12-12-73 
174 DIM 7(€3)sHC3)suC3) 
18% REM AVOGADRO'S NUMRER 
s | 194 LET V4=6e92R23 
24% REM ELECTRON REST ENERGY (MEV) 
91% LET Ra=.511 
22% REM CLASSICAL ELECTRON RADIUS ¢CM) 
239 LET D4=9.8E-13 
949 REM ALPHA PARTICLE CHARGE 
254 LET 74=2 . 
269 REM TARGET ATOM RADIUS (CM) 
279 LET R4=1F-8— 
ORF REM HARD SPHERE INTERACTION PRORARILITY 
29% LET Pas1F-5 
39% REM ATOMIC VUMRWR OF (€1=CUy 2=AGs 3=AU) 
314 LET 2¢19=29 
39% LET ZC HAZ 
334 LET 2¢€3)=79 
349 REM ATOMIC WEIGHT OF (1=CU, 2=AG> 3=AU) 
35% LET WC1)=63-5 
364A LET WC2)=197-9 
37% LET 8063)=197 
339 REM DENSITY OF (€1=CU> 2=AG> 3=AU) 
394 LET 4C€1)=8.96 
AVA LET HC2)=1%05 
419 LET H€3)=19.3 
hA=5e¢3 
S=3 


T 
x 
RT oT=1— 
4 
T 








SA=] 

= LA4A=19 

A779 PRINT’ INSTRUCTIONS C9SN0>5 1=YES) "5 

ASS INPUT I 

499 IF {=8 THEN 694 

594 IF I<>1 THEN 474 es 

514 PRINT'SCATRe PLOTS THREORETICAL PREDICTIONS OF THE FRACTION oF 
524 PRINTVALPHA PARTICLES SCATTERED AT ANGLES FROM 1a TO 184" 

5934 PRINT'DEGREES BY THREE DIFFERENT HEAVY ATOM MODELS." 

544 PRINTYINPUTS O&RE: ALPHA PARTICLE ENERGY CRETWEEN © AND 8 MEY)s’" 
994 PRINT'TARGET MATERIAL Cis=sCOPPER>, 29=SILVERs 3=GOLD), TARGET THICKNESS" 
969 PRINT'CRETWEEN «5 AND 2 MICROMETERS), DETECTOR SIZE" 

57% PRINT'CRETWEEN .5 AND 2 SQ CM)s AND DISTANCE FROM TARGET TOT 
939 PRINT'DETECTORS CBETWEEN 14 OND 49 CM)." | 











594 PRINT 

649 PRINT'ENERGY=""3 Ras 

619 INPUT I 

52% IF Cl-2)*(8-1)>=% THEN 659 

624 PRINT "OesPYERGY<=8" 

BAA GOTA 650 

654% LET FRa=1 : 

66% PRINT'MATERTAL="3 93 

6479 INPUT I 

BBA TR Cl=1)*C1-29*C1-3)=09 THEY 719 
69% PRINT =CUs P=AG, 324A" 

74% GOTO 66a 

710 LET S=!I 

729 PRINT "TYHICKNESS="5 Ts 

73% INPUT I 

TAM TF C105) C281) >=9 THEN 772 
754 PRINT". 5<=THICKVESS<=0" 

76% GOTO 72a 
779 LPP THT 
7343 LET T4s1E- 
79% PRINT'SIZEs" 
B99 INPUT I 

B14 LE C1-e5)# (9-1 )>=9 THEY Baa 
29% SPRINT". 5<=SI7texo" 

239 GOTO 79a 

24% LET Adal 

A549 PRINT'DISTANCK='"SLa43 

SKA INPUT T | 

RFY UF Cl=-19)*CAI-1)>=9 THEN 904 

R3% PRINT'LA<=DISTAYCRe=xa4a" 

R94 GOTO 854 

999 LET Last 

91% REM TARGET PARTICLES PER cmre 

932A LET N=T4*VU4eH¢(Sd7CS) 

934 REM DISTANCE OF CLOSEST APPROACH 

949 LET D=Z4*ZCS)*RA*DA/SEA 

954 REM MODEL CHOICE _ 

969 PRINT'MODEL Cl=HARD SPHERES» 2=THOMSON, 3=RUTHERFORD)"? 
979 INPUT M | 

Q9BA TF (M=$1)*£(M=$2)%(M-3)<>% THEY 969 

99% REM PREDICTION PLOTS 

1400 PRINT 

191% IF M>1] THEN 1459 

14926 PRINT"HARD SPHERE MODFI,"" 

1930 LET FsP4#N*R41O*AA/C ALL AIO) 

194% LET F1=8+19*(LOGCF)/LOG(1994+9)_ Ge 

1959 IF M<>2 THEN 1089 

1969 PRINT'THOMSON MODEL" 

1070 LET Fle? .. 

1989 IF M<3 THEN 1109 

1494 PRINT “RUTHERFORD MODEL" 

1194 PRINT E43;"MEV ALPHA'S “sTs "MICRON ‘ts. 

1119 IF S>1 THEN 1139 

1129 PRINT'"COPPFR': 

1130 IF S<>2 THEN 1159 

1149 PRINT'SILVERS 14 

1159 IF S<3 THEN 1173 | 


AAS 


4. 
i 
i. 
bid 

















SCATR2 (cont.) 


1144 
1179 
114A 
1199 
1209 
Leis 
1229 
1239 
1 2A 
1°56 
126423 
1279 
1 oR 
1290 
1 34e 
13194 
132A 
1333 
1349 
1354 
1363 
1374 
1389 
139% 
1A99 
141a 
1 4Aa0 
1439 
1444 
1450 
146% 
LATA 
1439 
149% 


PRINT" GOL D's. 
PRINT™ "5$A443"SQO CM COUNTERS "3L4a3"cm 
PRINT | 


AWAy" 


PRINT TABCI5)3 "FRACTION OF PARTICLES SCATTERED (LOG SCALE)" 


GOSUR 1314 

FOR 1214 TO 189 STEP 19 

LY M<a. TuEN 1279 

REM RUTHERFORD MODFIL 

LET Sl=SINC3.1416*1/360)1t4 
LET F=N*kD*D*kA4/C16*L4129%S1 ) 
Let P1=8+19*(LOGCF) /LOG(19)49) 
PRINT IS TARCTI$ Ts 

IF ABSCF1=389531 THEN 1379 
PRINT TARCFKL) 304" 

GOTO 1499 


PEN fe tm-9 1e-8 LE-7 lE=+§ 
PRINT. & LE-A LE= 3" 

PRINT “ANGLE: '*3 

PRINT TABC7>d3 7 * e eof * © eef 2 © eef 


PRIN Ge Fo. ® 
RETURN 

[F F1<7 THEN 12399 
PHINT TABC5733 
PRINT MCR ert ye 
NEXT 7 
GOSUS 
PRINT 
PRINT'NEXT PREDICTION CL=NEW MODEL, S=NEYW 
PRINT’ 3=END PROGRAM'S 

(INPUT 

IF X=1 THEN 964 

TF K=2 THEN 699 

IF K<>3 THEN 1439 

END 


ewe} * 


1344 


15 


ENERGY, 


LS" 5 





SCATR3 


REM SCATR3 - ALPHA PARTICLE TRAJECTORIES 
REM COPYRIGHT 1973 - STATE UNIVERSITY OF NEW YORK 

PM MODELS DEVELOPED SY Je RBRESMANs Ae CAGEIANOs AND De SCART, 
REM PROGRAMMED BY Ce LOSIKs Me WEISNERs De SCARL, AND Te 27.90" 
RM SCATER 1S 4 SET OF ASSOCIATED PROGRAMS CONTAINING SCATRI>» 
REM SCATR@, AND SCATR3 

REM LATEST REVISION: 12-12-73 

DIM 42€3)5RBC€6)9X%€6)24063915)9L69) 

rmeEM ELECTRON REST ENERGY CMEV) 

LET B4=.511 

REM VELOCITY OF LIGHT CFEMTOMETERS/SFEMTOSECOND) 

LET C4=s3E8 

REM CLASSTICSL ELECTRON RADIUS CFM) 


LET Dpa=2.28 


REM ALPHA PARTICLE CHARGE 
LET Z74=2 | 

REM ALPHA PARTICLE REST ENERGY CMEVD 

LET M4=3729 

REM TARGET ATOM RADIUS CFM) 

LET R4=2k5 

REM ATOMIC NUMRER OF (1=CU, 2=AG>, 3=4)) 

LET Z€1)=21 

LET ZC2)=A7 

LET 2€3)3=79 

LET F4=5e3 

LET S=3 

LET G=& 

FOR P=7 TO 9 

LET LCP y=4 

NEXT >? 

LET Nt=% 

PRINT 

PRINT 

PRINT"INSTRUCTIONS (S=NOs1=YESI"3 

INPUT I 

IF t=% THEN315 

IF I<>1 THEN 265 

PRINT'SCATR3 PLOTS THE TRAJECTORIES OF ALPHA" 
PRINT'PARTICLES SCATTERED FROM INDIVIDUAL YRAVY ATOMS." 


PRINT™INPUTS 4RE: ALPHA PARTICLE ENERGY CRBETWEEN 2 AND & MEV)," 


vege OS) ee MATERIAL Cl=COPPERs S=SILVERs 3=GOLD)," 

RINTVAND THREE DIFFERENT PERPENDICUL4SR DISTANCES RETWEEN THE" 
Shiga weioeNt ALPHA PATH AND THE CENTER OF THE TARGET ATOM." 
PRINT'ALPHA ENERGY="3 E43 
INPUT I 
LE 01-20% 08-19 Se9 7 
PRINT “@<=ALPHA ENE 
GOTO.315 
LET. BA4at 
PRINT'TSRGET MATERIALH=H"S S3 


BAF 


HEN 
eGye=sg" 
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SCATR3 (cont.) 


359 INPUT I 
355 IF €f-1)*(€3-1)>=9 THEN 379 
* 364 PRINT "“1=CUs 2=AG, 324" 
365 GOTO 345 
379 LET S=1I 
375 REM DELETE THE NEXT LINE TO MODIFY MAGNIFICATION EXPONENT 
38% GOTOA15 
385 PRINT™MAGNIFICATION=2-5E"3 G3" 9.5E"3 
399 INPUT I iS | 
395 IF (1-8)*C11-1)>=0 THEN 419 | 
AAS PRINT “&<=MAGe EXPONENT<=11" 
ASS GOTO 385 
ALA LETGsI 
415 LET G1i=141¢12-G) 
6 429 LET Ge=25«G1 
425 LET G3=3.2*6G1 
A389 LET G5=7*G63 
435 LET R4=R4/G63 
449 LET RA=R412 : 
ANS IF Ni<>3 THEN 465 
45% FOR P=1 TO 3. 
455 LET X¥CP=G1L*(2eP-298) 
46A NEXT P 
465 PRINT'CHOOSE X() RETWEEN'S -.8*G23 "AND" oSk&GOs "CP EMTOMETERS)" 
47% FOR P=1 TO 3 
475 PRINT ™XC'Ss Ps y="3xXCP)3 
48% INPUT 1 
485 IF .8*G2-ARSC(1I)<% THEN 475 
A9A LET XCPd=s] 
oT 





2 
, 


1, 
495 LET X(P+3)=1 
545 PRINT 
514 PRINT : 

515 REM RADIAL ACCELERATION CONSTANT 
528 LET 41574*705)*24*eD4eCAtosMg 
525 REM DISTANCE OF CLOSEST APPROACH 
534 LET D=Z4*Z7(S)*R4*kDA/EA 
535 REM INITIAL VELOCITY 
540 LET V4=C4*SOR(O*xE AMA) 
545 FOR P=1 TO 3 
SS FOR J=1 TO 15 
©... §$55 LET Atp..J=9 
564 NEXT J 
| 565 NEXT P 
6 579 FOR P=1TO03 

: 575 LF X€P)<>9 THEY 599 

5848 LET Q=3-14159 
585 GOTO 595 

59% LET O=ABSCO*XATNCD/C2*X(P)))) 

595 LET N1=SGNCXCP))*SGNC1.57979-0) 

699 FOR Y1=1TOS 

635 LET ACPsYLI=EINTCABSCXCP)/G1 +05) SGNCXCP) 426 
61% NEXT Y1 

615 ITF APSCN-1657979)>093 THENG39 

6249 LET ACPs Z=ACPs 8) /1994+25*SGNCXCP) 426 
625 GOTO 715 

639 IF Q<1.57979 THENG6G45 | 17 

&35 LET Y= | ; 





a, yt 








SCATR3 (cont.) 


649 GOTO 654 : 
G45 LET Y=ti5 : cI 
654 FOR Y1l=STOY STEP SGNC1.57979-0) 

655 LET X=XCPI+E3*CYV1-8)*4RSCSINCO)/C0S(Q))*V1 

569 TF Y1l=8 THEN 674 

665 LET ACPs YIIRACPSYLII/SIASF INTCARSCK/6G194+~.5)*SGNCX)4296 

679 IF ARSCX%)>G2 THEN 4685 

675 NEXT Yi 

6549 G0TO Fits | 

645 IFARSCG2-4RS(X)<=ARS(CEPQ-ABSCINTCACPs Y1-1))=-26)*G61) THENT7AS 

699% LET N=Y1-SGNC1-57%79-4) 

695 LET ACPsYLISCACPs Y1I-INTCACPs Y1929*%109 : 

744 GOTO 7123 a 
735 LET N=Yl 
TIA LET ACPesNISACPSNIAINTCACPs NN) 495% S56N6% 9426 

71S NEXT P 

720 PRINT TARCI3)3 "ALPHA ENERGY: "3 E43 "MEU" $s TARC 4303 "TARGET ''3 

725 PRINT CHARGES "3 2¢5S) 

739 PRINT 

T35 PHINTTARCI@I$S "XC LI="3 XC193 TARCIAIS PXCAI=""§ XCHI$ TARC HAS 9K ADH XC 2) 

T4% PRINT 

745 PRINT TASC 81S "CFEMTOMETERS)" 

794% PRINT TARCK)3"X="5 TARC1533-29*G1$ TARCO5)3 $1 AX GTS TARC 3572 atts 

755 PRINT TARC45)3 149*613 TARC55)324*«G!1 

4G9 PREV T'Y="5 

7165-G90SUB 935 


77% SOR Y= TG 15 = 
775 1% Y<>8 THEN 785 — 


405 LET Lez se6 
7385 IF ARSCY<-8)>R4 THEN 894 
194 LET LOSE INTC3e 2* SORCR3= CY¥-3)12)9405)4+296 


Vt CC8 S58) 
BAA ROR Pe} TO 3 
t LCP SIVTCACRs, ¥)) 
LOPtS = INTCCACEs YIKLOCP))*1494+.5) 


ATS NEST PE 

G4 PRINT GA¥CY=S)3 TARCIBIS"INS 
895 LETI=62 

B34 FOR Pal Ta: 9 


S59 NEXT P 

855 IP J=62THENI 53 

$63 LET L€I)=9 

865 PRINT TA8C14+J3)3 

AIA IF Y¥<>S8THEN $89 

S75 IF J=26 THEN 9123 

S84 IF - €1l$1)*C1=4)<58 THEN 899 
‘heres Cet ey METS 

BIS TF Cl-+2)*(1-5)<sa THEN 999 
SIO CETL ERs 


a ee 
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SCATR3 (cont. ) 


Jag 
QAS 
914 
915 
929 
95 
93% 
935 
9AG 
P45 
95a 
955 
959 
965 
979 
975 
998A 
985 
999 
995 
1499 
1ZAaS 
1419 
1415 
142% 
1AL5 
1 A308 
1435 
1440 
14A5 
145% 


TRE €1-3)*¢CI[-6)<>9 
PHI 33 


THEN 916 


“TF [<>7 THEN 926 


PRINT "N''s 
IF I1<8 THEN 8285 
PRINT "Rt; 
GOTO 825° 
PRINT 
PRINT 
RETURN 
PRINT 
IFY=15 THEN 965 
PRINT 
NEXT Y¥ 
GOSUR 935 
FOSP=1TO3 
PRINTTARC232*CPH1dIISMANGLEC'S P3 ty="'3 
TF X(P)<>9 THEN 10249 
PRINT 1883 
GOTO1395 
PRINT INTCARS( 369/36 1LALSOXATINOD/SCE*XCPIIII“ELABDIS1AAS 
YEXTP 
PRINT 
PRINT 
PRINT™L=NEY 
INPUT I 
PPINT 
IF I=1 THEN 473 
IF f=2 THEN 315 
IF I<>3 THEY 1929 
END 


—=— se =e — ee — ee on ee oe <a 
we 


DISTANCES» 2=NEW ENERGY OR TARGET, 
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3=END PROGRAM'S 


IV. SCATER NOMENCLATURE 











So a Avogadro's number = 6.02 E23 atoms/mole 
B4 = Electron rest energy = .511 MeV 
C4 = Velocity of light = 3 Eo meters/second = 3 E8 femto- 
. meters /femtosecond 
D4 = Classical electron radius = 2.8 E-13cm = 2.8 femtometers 
Z(1) = Atomic number of copper = 29 ' 
Z(2) = Atomic number of silver = 47 . 
Z(3) = Atomic number of gold = 79 
W(1) = Atomic weight of copper = 63.5 amu 
W(2) = Atomic weight of silver = 107.9 amu 
W(3) = Atomic weight of gold = 197 amu 
H(1) = Density of copper = 8.96 g/cm? 
H(2) = Density of silver = 10.5 g/cm? 
H(3) = Density of gold = 19,3 g/cm 
E4 = Alpha particle energy (MeV) 
Z4 = Alpha particle charge = 2 @ 
M4 = Alpha particle rest energy = 3720 MeV ae 
B = Number of beam particles 
Bl = Log of number of beam particles 
R4 = Target atom radius = 2 E-8cm = 200,000 femtometers 
T4 = Film thickness (cm) 
t = Film thickness(micrometers) 
P4 = Hard-sphere interaction probability = 1 E-5 
A4 = Counter area (cm) 
L4 = Counter distance from target (cm) 
M = Model index: l=spheres; 2=Thomson; 3=Rutherford 
S = Material index: l=copper; 2=silver; 3=gold 
D = Distance of closest approach : 
N = Number of target particles per square cm. 





ner 
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